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Abstract 

We compare the probability distributions and Binder cumulants of the overlap in 
the 3D Ising spin glass with those of the magnetization in the ferromagnetic 2D XY 
model. We analyze similarities and differences. Evidence for the existence of a phase 
transition in the spin glass model is obtained thanks to the crossing of the Binder 
cumulant. We show that the behavior of the XY model is fully compatible with the 
Kosterlitz-Thouless scenario. Finite size effects have to be dealt with by using great 
care in order to discern among two very different physical pictures that can look very 
similar if analyzed without large attention. 



1 Introduction 



The issue of the existence of a phase transition in the three dimensional Ising spin glass has 
been a hard and difficult problem for more than two decades (see for example |jl|, and 
references therein). Today there are clear numerical evidences favoring the existence of a 
low temperature broken phase 0, ^, |], 0], but a deeper understanding of the underlying 
physics still lacks. It is clear, for example, that one is very close to the lower critical 
dimension (LCD), but understanding the details of the influence of such effect is highly 
non-trivial. 

Determining for example the infinite volume limit of the Edward- Anderson order pa- 
rameter (gEA) has been beyond reach till very recently, and the existence of the phase 
transition (both in 3 and 4 dimensions) has been established by exhibiting the crossing of 
the finite size Binder parameter. One was able to show (for the 3D case see P, ^) that 
curves of QhiT), g^^iiT) as a function of T would cross at T]-^), but it was impossible to 
determine the non-trivial limit of giiT) for L ^ oo at T < Tc (and in the same way it 
was impossible to determine the large volume limit of ^ea)- Only in the most recent pe- 
riod off-equilibrium techniques [0 and equilibrium simulations based on parallel tempering 
P] have allowed a statistically significant determination of the 4D infinite volume order 
parameter gEA- 

We start here by noticing that the behavior of the Binder parameter giici) in the 3D 
spin glass is very reminiscent of the one one finds in the 2D XY model (without quenched 
disorder), ^'^(m) (here m is the magnetization). Even for quite large lattices the curves for 
different lattice volumes are well split in the high T phase, but seem to merge better than 
cross at low T. Only on very large lattices one can exhibit a non-ambiguous (but always 
very small) crossing |^. The XY 2D model shows that the order parameter is zero in 
the thermodynamical limit only very slowly when increasing L. 

The same kind of effect could be appearing in the 3D Edwards-Anderson spin glass, 
and in order to be sure one is dealing with a real phase transition with a non-zero order 
parameter one has to be very careful, and to show that is keeping under control possible 
contamination. That is why we have decided to run a comprehensive comparison of the 
order parameter distributions for the 3D Edwards- Anderson spin glass and for the 2D XY 
model. A detailed paper by Binder 0, containing a study of the distribution functions 
for the Ising model, can be considered a methodological prototype to this kind of analysis, 
and can be used as a nice introduction to the finite size scaling techniques and ideas used 
in this setting. 

Let us start by reminding the reader about some main points concerning the definition 
of the lower critical dimension. The lower and the upper critical dimensions [di and du 
respectively) are important in qualifying a statistical system, du is the minimal dimension 
where mean field predictions hold (apart from logarithmic corrections), while the LCD, 
di, is the maximal dimension where the finite T phase transition disappears. A typical 
example is the usual Ising model, with di = 1 and (i„ = 4 |Tr!|. 



Since a (p^ term appears in the effective Hamiltonian of spin glasses (see [Tl|] and 
references therein) one expects that the upper critical dimension is du = 6. One of the 
possible ways to determine the lower critical dimension is based on the determination of 
the critical exponent r]. One starts from the two points correlation function at the critical 
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point, T = Tc, that for \x — y\ — oo behaves as 

(0(f)0(f))^if-yr('^-2+'') . (1) 

The LCD is defined by 

di-2 + 7]{di) = , (2) 
i.e. by the fact that there is no power law decay of the two point correlation function at 



the (T = 0) critical point. A (replica-symmetric) e-expansion computation gives 

7] = --e+ 1.2593e2 + 2.5367e3 , (3) 
3 

where e = (6 — c?). At order e one is getting the promising estimate di = 3, that collapses 
when including the higher order contributions, that do not allow any real solution for 
di < 6. It is clear that because of one of the many reasons that could give troubles (for 
example replica symmetric breaking and poor convergence of the e-expansion) here the 
e-expansion is not helping in determining the LCD. 

Equation (^ allows an estimate of di based on numerical estimates of the rj exponent. 
In four dimensions, with Gaussian couplings, one finds rj = —0.35 ± 0.05 0, while in 3D 
rj = —0.40 ± 0.05 [0. The variation of t] with d is small, and it seems safe to estimate 
di ~ 2.5. Even if this result is somehow peculiar (since in the field theoretical approach 
|TT|| one does not see any trace of propagators with non-integer powers) it is confirmed by 



a mean field based analysis |]13[ where one builds up an interface and looks at its behavior. 
This mean-field computation gives di = 2.5, in excellent agreement with the numerical 
estimate. 

Numerical simulations in 3D |P, ^, |^ have now shown clearly that there is a finite T 
phase transition, i.e. that di < 3. The broken phase is mean field like, and understanding 
more details about it will be the goal of this paper. It is also well established that in 2D 
one finds a T = phase transition (see and references therein). Summarizing, from 
state of the art numerical simulations one can deduce that 2 < di < 3. 

Also the fact that d^ = Q is well supported by numerical results ITl, IT3. In 6d one 



determines with good accuracy mean field exponents (7 = 1, /5 = 1 and z = A), with 
logarithmic corrections (that have been detected in the equilibrium simulations). 

Here we will try to shed more light on the difficult numerical simulations of the 3D 
Edwards Anderson spin glass. The main problem is probably in the fact that the system 
is very close to its LCD. So the apparent merging of the Binder parameters in the low 
T region, that has only recently been disentangled to show a significant crossing [0, is 
dramatically reminiscent of the one one can observe in the case of a Kosterlitz-Thouless 
transition. We will try here to learn more about the effects of an anomalous situation like 
the Kosterlitz Thouless (KT) one, by looking for example to the Binder cumulant and to 
the overlap probability distribution P{q). To do this we will discuss in same detail the 
structure of the order parameter probability distribution in the 2D XY model without 
disorder. We will stress how similar to the 3D spin glass things are at a first level of 
analysis, and where the relevant differences can be found. It is also remarkable that the 



pure 2D XY model has a peculiar aging behavior [16|: aging is one of the crucial features 



of spin glass systems, and its qualification is of large importance. 
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In the next section we will define our models, the physical observable quantities, and we 
will give details about our numerical simulations. In section we discuss our results, by 
following in parallel the 2D XY model without disorder and the ?>D spin glass. In section 
(|) we draw our conclusions. 



2 Models, Observables and Simulations 

We have studied the two dimensional XY model on a squared lattice. The volume is 
denoted hy V = L"^, the hamiltonian is 

n = - cos(0. - <Py) , (4) 

<x,y> 

where < x,y > denotes a sum over nearest neighbor site pairs, is a continuous real 
variable, and periodic boundary conditions are imposed on the system. This model shows 
an infinite orc?er phase transition (with (3^ ~ l-H) [|T3' Kosterlitz-Thouless transition. 



In according with the Mermin- Wagner theorem there cannot be non-zero order pa- 



rameters: the magnetization in the thermodynamical limit is zero for all T > 0. The KT 
transition is characterized by a change in the behavior of the two point correlation func- 
tion, which goes from the exponential decay of the high temperature phase to the algebraic 
decay of the low temperature phase. The whole low temperature phase (/? > Pc) is critical 
(the correlation length is infinite). 

To simulate this model we have used the Wolff single cluster algorithm |TB[. The 
simulations have been run at five different values of /3 in the low temperature phase: 
P = 1.3, 1.4, 1.5, 1.7, 2.0. For each value of P we have used the lattice sizes L = 8, 16, 32, 
64, 128, 256. For each value of (/?, L) we have used 200, 000 iterations of the single cluster 
algorithm, discarding the first half for thermalization. The total CPU time required has 
been approximately one month on a 100 MHz Pentium based computer. 

We have measured the probability distributions of 



1 

m^ = — 
V 



Re ^exp(i0^) 

X 



(5) 



that we denote as Pi (mi). We call m™'^^ the value of mi where Pi (mi) is maximum and 
takes the value P™^^ = Max[Pi(mi)]. We have looked in detail to the first and second 
moments of Pi (mi), (mi) and {mf). We have also computed the Binder cumulant of the 
Pi (mi) distribution: 

- m . («) 

{miyj 

At low T one can study the XY model by using the spin wave approximation, that neglects 
the role of vortices (since they are suppressed at low T). In the T limit all the spins 
point in the same direction, and {6) is uniformly distributed, so that 

(m?) = — r dOcoi^^e , (7) 




2tx Jo 
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and the Binder cumulant at T = has the value 

B,{T = ^) = \. (8) 

To determine the relevant scaling behavior we use the fact that, for the XY model, 
X — L^~^*^"^) and {m?) = x/L"^-, where in the spin wave approximation 

'KT) = I , (9) 

is the anomalous dimension of the field. Moreover, since Pi (mi) is a probability distribu- 
tion, normalized to one, with non zero maximum value (m™^^) (at least for finite values of 
the lattice sizes) and with (mi) ~ L^"^ (at finite temperatures) we have thatQ 

independently of the lattice size, L. Since 

m^" ^ (mi) ^ (m2)i/2 , 

we conclude that 



mf*^^ 


^ ^-r?(T)/2 


(mi) 




{ml) 


^ T ~V(T) 
— ^ ) 


pmax 
1 


- /,'7(T')/2 _ 



(11) 

The other model we have studied is the three-dimensional Ising spin glass with quenched 
random couplings J distributed with a Gaussian law. The Hamiltonian is 

71 = - ^^■k,^3 > (12) 

<i,3> 

where the spin are defined on a three-dimensional cubic lattice and < i,j > denotes a sum 
over nearest neighbor pairs. 

As usual we have simulated two real replicas (cr and r) with the same quenched 
couplings, and we have measured the overlap 

q{cr,T) = ^Y.^i'^i ^ (13) 

i 

and its probability distribution 



P(g) = (5(g-g((7,r)) (14) 

where, as usual, we denote the thermal average with ((■ ■ ■)), and the average over the 
disorder distribution with (■■■). The Binder cumulant of the probability distribution P{q) 
is 



^In the rest of the paper the symbol A o:^ B means A = 0{B). 
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9 = 1; 



(15) 



We have run L = 4,6,8,10,12 and 16 lattices with 2048, 2560, 512, 512, 2048 and 500 
samples respectively. We have used the supercomputer APE-100 
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For simulating the spin glass model we have used the simple tempering method for 
small lattices {L < 10), and the parallel tempering scheme for large lattices (L > 12) (see 
|pO| , ^ 1^ and references therein). Thanks to that we have kept under control the level 
of thermalization reached by the system, that is in all cases very good (for a discussion of 

PH). We have checked that the equalities estabhshed 
] hold for our results, and that the P{q) is well 



the standard criteria of control see 
numerically in M , and proven by Guerra 



symmetric, supporting the reach of full thermalization. 



3 Results 




Figure 1: Binder cumulant for the 3D Ising spin glass. On the right, from top to bottom, 
curves and data points are for L = 4, 6, 8, 10, 12 and 16. 

In figures (0) and @ we show the Binder cumulant and the probability distribution of 
the 3D Ising spin glass. 

Let us discuss first the Binder cumulant. In figure (|l|) there are two different regions. 
In a high temperature region curves corresponding to different lattice sizes are clearly split 
(they tend to zero in the thermodynamical limit). At small temperature (i.e. for f3 larger 
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than f3^^ ^ 1-0), on small lattice sizes (up to L = 10) curves coalesce, within our small 
error bars, in one. It is interesting to notice that defining a Binder cumulant based on 
three different replicas pSf allows a somehow easier determination of the critical behavior. 




Figure 2: Probability distribution of the overlap, P{q), for the 3-D Ising spin glass. T = 0.7 
From right to left curves and data points are for L = 4,6,8,10 and 16. 

Only when thermalizing a L = 16 lattice (quite large for current standards, and im- 
possible to thermalize deep in the critical region without the use of parallel tempering 
0) one is able to exhibit a clear crossing between, for example, the L = 8 curve and the 
L = 16 curve. This implies the existence of a phase transition at finite temperature with 
a non-zero order parameter, gEA 7^ (see Kawashima and Young for the model with 
quenched binary couplings, J = ±1). 

We can compare figure (Q) with figure (^, where we show our numerical results for the 
Binder cumulant, Bi, for the 2D XY model. Up to L = 10 the XY model and the 3D Ising 
spin glass have a very similar behavior: again, within error bars, in the low temperature 
region all the curves for different lattice sizes collapse in a single curve, without any visible 
sign of finite size effects. 

The behavior of the full probability distribution of the order parameter (mi for the 2D 
XY model and q for the 3D Ising spin glass) is very similar. Figures (§) and (^, where 
we show Pi (mi) ai j3 = 1.3 and (3 = 2.0 respectively, can be compared to the analogous 
figure for the spin glass P{q), (@). The overall shapes are very similar. The peak shifts to 
the left, and in both cases -P(O) looks constant in our statistical precision. 

We give in table (P, at /? = 1.3, the expectation values of the observables shown figures 
(|) and (^). By fitting these values by using a single power fit we find 
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P 

Figure 3: Binder cumulant, Bi, for the 2D XY model. Prom top to bottom, L = 4, 8 and 
16. 



3 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
+ L=8 


1 1 1 1 1 1 










° L=64 f /[ jfe 




■ '-=128 / JfVn 




V L=256 // /I f 1/ 








1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ivA 





0.2 0.4 0.6 0.8 1.0 



Figure 4: Probability distribution for the 2D XY model, Pi (mi), ai P — 1.3. 
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Figure 5: Probability distribution for the 2D XY model, Pi (mi), at /3 = 2.0. 



^max ^ ^-(0.08±0.01) 
(mi) ~ /,-(0.09±0.01) ^ 

{m\) ~ L~^o.i7±om) ^ 

pmax ^ ^+{0.m±Om) j^j^g^ 

These results are in remarkable agreement among them and in good agreement with the 
spin wave exact value rj{[3 = 1.3) = 0.12. Corrections due to vortices are equivalent to a 
higher effective temperature that undergoes here a 30% shift. 

We have also established that a power fit to a non-zero infinite volume order parameter 
of the form 

^max(^) = m^^"(00) + (17) 

with m™'^^(oo) different from zero and A and B constant is excluded by the data. 

Figures (^ and (|^) are interesting: they show a finite size non-trivial behavior that 
we know, from theoretical ideas (the Mermin- Wagner theorem), and from the analysis of 
the numerical data, will converge to a zero centered delta function limiting probability 
distribution in the infinite volume limit. This is the point we want to stress. Since the 3D 
spin glass has a very similar behavior (and even for the Ad model, where the crossing of 
the Binder cumulant is clear, it is non-trivial to show that qea tends to a non-zero limit) 
it is crucial to understand where differences are. 
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We also want to stress that Pi (mi) shows a clear plateau, roughly L-independent, close 
to the mi ~ region. The plateau height grows with the lattice size (in a statistically 
significant way in our numerical data, see (|I|)). This is one of the interesting results of 
this note: the infinite volume 2D XY mi delta function is constructed from the increasing 
finite volumes by a finite mi peak that shifts towards mi ^ 0, and by a plateau in the 
mi ~ region that slowly increases with the lattice size, to eventually match the peak in 
the mi = delta function. 

We have repeated this analysis for the overlap probability distribution P{q) of the 
?)D spin glass (see figure (0)). The best scaling fit of the peak position, g^, where the 
probability distribution is maximum, by a power law (the data are in table @)) gives 

qj^ = (0.70 ± 0.02) + (1.6 ± 0.7)L-(^-^±°-^) , (18) 

where T = 0.7. In this fit we have used all lattice volumes (L < 16). The fit had a 
X^/dof = 0.15. This thermodynamical value we get for ^ea is close to the value that has 
been extracted from an off-equilibrium simulation (g ~ 0.7) 0. The best (two parameter) 
fit obtained by fixing qm = 0.7 (considered as an input from the dynamical simulations) 
gives compatible results with smaller errors. In figure (P) we show the qm data versus L~^'^ 
(see also table (^), and the curve from the best (two parameter) fit. 

From the numerical data for the 3D spin glass (that are from a state of the art large 
scale numerical simulation) we cannot exclude the possibility of qea = in the infinite 
volume limit. We find that the best fit (that uses in this case only L > Q data) 

gM = (1.0 ±0.1)L-(°-i2±o.o2) (19) 

is very good. So, even if the scenario of a non-zero overlap is favored (the static value is 
equal to the dynamic one, the exponent of a decay to g = is very small) in the 3D case we 
cannot use this limit to be sure of the existence of a phase transition with a non-zero order 
parameter (in Ad recent high statistics data allow to establish this evidence ||^). The safe 
evidence for the existence of a phase transition in the 3D spin glass relies in this moment 
on the statistically significant crossing of the finite L Binder cumulant |^, that makes 
visible fine details of the equilibrium probability distribution. 

Also the behavior of -P(O) turns out to be the potential source of many ambiguities. 
We have seen that in the XY model it grows very slowly with the lattice size, in order to 
asymptotically contribute to the mi = delta function. Also in the 3D spin glass, where 
in a mean-field like broken phase we expect a finite limit for P(0) we observe a constant 
plateau with a (non necessarily statistically significant) growth for L = 16. This behavior 
contributes to falsify the droplet model picture, where one would expect -P(O) to decrease 
with the lattice size. 

4 Conclusions 

We have shown that the two dimensional ferromagnetic XY model and the three dimen- 
sional Ising spin glass finite volume order parameter probability distributions behave very 
similarly. The Binder cumulants on small lattice volumes show a similar merging at low T. 
Only on large lattices the 3D spin glass exhibits a crossing typical of a phase transition. 
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Figure 6: Value of the overlap Qm such that P{q) is maximum {ZD Ising spin glass). The 
continuous line is the fit described in the text. T = 0.7. 
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(mi) 


{ml) 




A(o) 


16 


0.74(1) 


0.485(6) 


0.283(4) 


2.64(9) 


0.83(5) 


32 


0.70(1) 


0.456(5) 


0.260(3) 


2.79(8) 


0.93(2) 


64 


0.66(1) 


0.431(5) 


0.231(3) 


3.00(9) 


0.97(6) 


128 


0.62(1) 


0.405(4) 


0.205(2) 


3.08(7) 


0.94(9) 


256 


0.60(1) 


0.382(5) 


0.183(3) 


3.36(11) 


1.09(6) 



Table 1: Numerical data for the 2D XY model, (3 — 1.3. See the text for more details. 



L 


qu 


P(0) 


4 


0.91(1) 


0.398(3) 


6 


0.81(1) 


0.376(5) 


8 


0.77(1) 


0.39(2) 


10 


0.76(1) 


0.39(2) 


16 


0.72(1) 


0.49(7) 



Table 2: Numerical data for the 3D Ising spin glass, T — 0.7. See the text for more details. 
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The results that we have discussed for the XY model are completely compatible with the 
KT predictions. We have analyzed the finite volume behavior of the peak of the finite vol- 
ume order parameter probabihty distribution. In the case of the XY model the preferred 
limit is zero. In the spin glass case the preferred value is non-zero, and compatible with 
an off-equilibrium estimate, but from the present data one cannot rule out the possibility 
of the position of the peak going to zero in the infinite volume limit. 

We have also established that P{0) in the KT scenario has a finite volume non-zero 
value, that increases in the infinite volume limit. In finite volume one can then exhibit 
probability distributions with the same shape of that of a finite dimensional spin glass that 
have a trivial thermodynamic limit (a delta function in the origin). Analyzing finite size 
effects is crucial before reaching conclusions about the critical behavior. In the 3D spin 
glass good evidence for the existence of a mean field like phase transition is based on a 
dynamical determination of the Edward- Anderson order parameter and on the crossing of 
the Binder cumulant on large lattices, but determining with good precision the shape of 
P{q) on large lattice sizes will be important for making more details of the critical behavior 
crystal clear. 
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